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I. INTRODUCTION 



Because of the non-perturbative nature of the problem and the complexities of the re- 
action mechanisms, the calculations of meson-baryon reactions within the framework of 
relativistic quantum field theory can not be done exactly. Thus the progress we can make 
now and in the forseeable future is to construct manageable reaction models for analyzing the 
data. Each model involves some approximations and assumptions. For understanding the 
information extracted from the data, such as the electromagnetic form factors of the N-N* 
transitions, it is necessary to examine this phenomenological aspect of the employed models. 
This is the objective of this paper concerning the approaches based on (1) three-dimensional 
reductions (2) method of unitary transformation [2I, [s], 0], and (3) time-ordered per- 
turbation theory 0| . For simplicity, we will only consider the case of single-channel irN 
scattering. This is sufficient to reveal the differences between these approaches. We will also 
examine the relations between these approaches and the approach based on the dispersion 
relations of the S-matrix theory. As is well documented 0], these two theoretical frameworks 
are based on very different theoretical considerations and hence there is no compelling the- 
oretical reasons to favor one of them in developing phenomenological models to analyze the 
data. 

The considered three methods for constructing ttN models are explained and analyzed in 
sections II-IV. In each section we will also address their relations with the approach based 
on dispersion relations. In section 5, we give a summary. 

II. THREE-DIMENSIONAL REDUCTIONS 

To illustrate the derivations of three-dimensional equations for ttN scattering from rela- 
tivistic quantum field theory, it is sufficient to consider a simple ttNN interaction Lagrangian 
density 



where tplx) and denote respectively the nucleon and pion fields and Fq is a bare nNN 
vertex, such as Fq = ig'j^ in the familiar pseudo-scalar coupling. By using the standard 
method [8|, it is straightforward to derive from Eq. ([T]) the Bethe-Salpeter equation for 
ttN scattering and the one-nucleon propagator. In momentum space, the resulting Bethe- 
Salpeter equation can be written as 



where k and P are respectively the relative and total momenta defined by the nucleon 
momentum p and pion momentum q 



Here f]N{y) and T].,^{y) can be any function of a chosen parameter y with the condition 



(1) 




(2) 



P = p + q, 

k = vAy)p~ VNiy)q- 



V7r{y) + VN{y) = 1 • 



(3) 
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Obviously we have from the above definitions that 



P = VN{y)P + k, 

q = r]^{y)P-k. (4) 

In analogy to the nonrelativistic form, they are often chosen as: rjj^j = mj'q/{mT^ + mAr) and 
r/^ = m.,r/{^TT + '^Af)- The choice of the rj's is irrelevant to the derivation presented below 
in this section provided that Eq. Q is satisfied. 

Note that T in Eq. ([2]) is the "amputated" invariant amplitude and is related to the 
7tN S-matrix by S* oc uTu with u denoting the nucleon spinor. The driving term B in 
Eq. ([2]) is the sum of all two-particle irreducible amplitudes, and G is the product of the 
pion propagator D.„{q) and the nucleon propagator Sn{p)- In the low energy region, we 
neglect the dressing of pion propagator and simply set 

DA(i) = - — (5) 

where is the physical pion mass. The nucleon propagator can be written as 

^^^^ i;^ - m% - T.j^^i^p'^) + ie' ^ ^ 

where m% is the bare nucleon mass and the nucleon self energy operator Sat is defined by 

Siv(p') = / d''kToG{k;p)t{k-p). (7) 

The dressed vertex function F on the right hand side of Eq. ([7]) depends on the vrA^ Bethe- 
Salpeter amplitude 

r{k; P) = To + J d^k' To G{k'; P) T{k', k; P) . (8) 

It is only possible in practice to consider the leading term of B of Eq. ([2]). For the simple 
Lagrangian Eq. ([1]) the leading term consists of the direct and crossed diagrams 

B{k, k'] P) = B'^"^ {k, k'; P) + (A;, k'; P) , (9) 

where 

B^'^\k,k'-P) = Fo . ^ \ , . Fo, (10) 

B'^\k,k'-P) = To . p \ ^ . To, (11) 

with P = [r]N{y) - r]M]P + k + k'. 

Equations (|2I)- (ITT1) form a closed set of coupled equations for determining the dressed 
nucleon propagator of Eq. (jS]) and the vrA^ Bethe-Salpeter amplitude of Eq. ([2]). It is impor- 
tant to note here that this is a drastic simplification of the original field theoretical problem 
defined by the Lagrangian ([1]). However, it is still very difficult to solve this highly nonlinear 
problem exactly. For practical applications further approximations are usually introduced. 
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The first step is to define the physical nucleoli mass by imposing the condition that the 
dressed nucleon propagator should have the limit 



SNip)\p2^m% -7 , (12) 

where m^v is the physical nucleon mass. This means that the self-energy in the nucleon 
propagator Eq. is constrained by the condition 

m% + EN{m%) = rriN. (13) 

The next step is to assume that the p-dependence of the nucleon self-energy is weak and we 
can use the condition Eq. (fT3l) to set m% + T,{p'^) ~ m% + 'E{m'lf) = rriN. This approximation 
greatly simplifies the nonlinearity of the problem, since the full t^N propagator G in Eqs. 
(2), (7) and (8) then takes the following simple form 

G{k-P) = - V--. (14) 

The next commonly used approximation is to reduce the dimensionality of the above 
integral equations from four to three. The resulting models will be called 3dBS models in 
the following discussions. There exists extensive literature on this subject, as reviewed in 
Ref. [if. A rather complete numerical study of several of these 3dBS models was given 
in Ref. We therefore will not get into these subjects. Instead we will focus on how these 
models are interpreted in the actual analysis of nN data. 

Let us specifically consider a 3dBS model derived from using the three-dimensional re- 
duction method of Kadyshevsky [13]. In the cm. system, P = {E,6), we first define 



v{k\k,E) = -^./-^^ , ^ UrB{k',k,E)urJ-^^ , ^ = (16) 



where the time components of the momenta of the BS amplitude T and interaction B in 
Eq. (I2| have been fixed by the chosen three-dimensional reduction, as explicitly given in 
Ref. (9|. The resulting three-dimensional scattering equation can then be cast into the 
following form 

t{k',k,E) =v{k',k,E) + J dk"v{k',k",E)g{k",E)t{k",k,E) (17) 
where the propagator is 

'^''^''^ = E-E,ik)-EAk)+^e■ ^'^^ 

For simplicity, we only consider the case that the interaction v is derived from s- and u- 
channel mechanisms defined by Eqs.fl9|)- f|TT|) . The difference between 3dBS models whose 
scattering equation can be cast into the form of Eq. f|T71) is in the expression of the driving 
term v in Eq. (fT7|l . It is important to note that the only common condition these 3dBS 
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models satisfy exactly is the unitarity condition, as explained in Ref. [jj. Thus it is not 
surprising that their amplitudes have different analytic structure in complex-E plane. 
To simplify the presentation, we use the following operator form for Eq. ([1] 



t{E) = v{E) + v{E)g{E)t{E) (19) 
For most of the 3dBS models, one gets the following form of v{E) in the Pu partial wave 

v{E) = vP°^^{E)+v''\E) (20) 

with 

-'"''(E) = rl-r^To (21) 

where m% is the bare mass of the starting Lagrangian. The form of the second term v^^{E) 
of Eg. (1201) is not important in the the following discussions. 

Substituting Eqs. (I2U|) - (I^ into Eq. ([T^ . we then get the following well known form 

t{E) = t'\E) + — rt(ij;)r(E) 

where 

t^3{E) = v^3{E)+v^3{E)g{E)t^3{E), 



T{E) = To[l + g{E)t'^{E)], 

En{E) = {No\Tog{E)r\E)\No). (23) 

In the above equations \Nq) is the bare one-nucleon state, "Ej^ilE) and F are the three- 
dimensional forms of Eqs ([7]), ([H]), respectively. The condition Eq. (IT^ then becomes 

m% + ^^{171^) = rriN (24) 

which defines the bare mass from the physical mass rriN- It is common [ll] to use 
Eq. (l22l) at the E — > m^v position to define the dressed vertex T{E) in terms of physical 
coupling constant. To get this relation, we first note that at E rriN, the self-energy T^^^E) 
can be expanded as 

En{E) = S^(mjv) + {E- mjv)Si(m^) + • • • (25) 

where 



dE 

By using the above relations, Eq. (122|) becomes 



(26) 



tiE)\E^^^ = t'n^iv) + [Zl^'r^mN)]- [Zl^'VimN)] (27) 

£/ — TUn 
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with 

Z2-i = l-Si(m^). (28) 

1 /2 

The renormahzed vertex [Zg T{m^)] is then used to define the physical couphng. If the bare 

r(0) _ 

vertex Fq is chosen to be the usual pseudo-vector coupling Lj = ip'j^'jf^ipd^cl), the above 

procedures relate the renormahzed coupling constant fnNN to the bare coupling constant 
AO) 

J-kNNi 

Unn = + g{mN)t'^{mN)]Zl'\ (29) 

The renormahzed coupling constant is identified with the empirical value g^j^j^/Air = 
(2m^/m02(^V/47r) = 14.3. 

Here we point out that the above procedure implies an interpretation where the physical 
nucleon is made of a bare core Nq and pion cloud. To illustrate this, it is sufficient to 
consider the case when ^('•s) = 0. If we set {No\To\k) = To{k), Eq. ([22]) has the following 
analytic form, 



where 



t^k\k,E)=rUk) ^_^^,,_^,^^^ T„(k) (30) 



The nucleon pole condition Eq.(24) then becomes 

= + S^(mAr) (33) 

{k'\t{E)\k)\E^m. = m') ^ fo(fc) (34) 

-C/ — 777, TV 

where the renormahzed vertex function is 

fo(A;) = Z'o^^Toik) (35) 

with 

Zq-i = 1-S?(m;v) (36) 

= ^ + /^^^' (m^-K(l)-E^(A:))^ • (^^^ 

It is interesting to note that the above nucleon pole conditions can be reproduced by assum- 
ing that the structure of the nucleon can be described by the following mass operator h in 
a subspace spanned by the state \No) and \k) for the vrA^ state 

h = ho + To (38) 
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with 



ho\No) = m%\No) (39) 

ho\k) = {E^{k) + E^{k))\k) . (40) 

We assume that the physical nucleoli state is defined by 

h\N) = m^lN), (41) 



\N) = z, 



1/2 




A^o) + / edkf{k)\k) 



(42) 



The basis states are normalized as (iVolA'o) = 1, {k\k') = k'^6{k — k'). Projecting Eq. (14T|) 
from the left onto |A^o) and \k), we then get 

mo + J k'^dkf{k)To{k) = tun (43) 

[E^ik) + E^{k)]f{k) + T;{k) = m^fik) (44) 
The normalization condition (A^|A^) = 1 gives 

zo = l + J k^dk\f{k)\^ (45) 

From Eq. flH|) we have the solution 

^^^^ ^ mr, - Er^ik) - EM ' ^^^^ 

Substituting Eq. fH6|) into Eq. fj43l) . we then get exactly the nucleon pole condition Eq. fl33l) . 
Within the model defined by the mass operator Eq. pTl) . the physical ttNN vertex can be 
calculated from using Eq. fH2|) 

r^(fc) = {N\To\k) = z'J'Voik) (47) 

Substituting Eq. ( H6|) into Eq. ( 145|) . we find that 2o is exactly equal to Zq of Eq. ( 1371) and 
hence Tj^{k) is exactly the renormalized vertex Tq of Eq. (1551) . 

The above simple model illustrated that the nucleon pole conditions Eqs (33)-(37) can 
be related to the substructure of the nucleon. If we write Eq. (H2l) as an operator form 
|A^) = lA^o) + f\7^N) and iterate it, we then get 

\N) = zl^'[\No) + /IvriVo) + //IvrvriVo) + ///IvrvrvriVo) + ■ ■ ■ . (48) 

This illustrates that the usual procedure of requiring the vriV amplitudes to have a nucleon 
pole implies that the physical nucleon is made of a A'^o core and meson cloud. 

To be consistent, one in principle should also replace E^ik) in the propagator Eq. (ITS]) 
by an expression which is related to bare mass m^^r and the matrix element of the self energy 
S^, defined in Eq. (|3T|) . But this complicates the unitarity condition for the scattering 
amplitude t defined by Eq. (|T71) . To do it properly, one needs to also consider the rniN 
unitarity condition since the self energy contains ■nN intermediate state. 
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Let us stress that the above observation has a connection with the vrA^ scattering equation 
derived by Aaron and Amado and Young [l^ (AAY) using the three-dimensional reduction 



of Blankenbecler and Sugar [13j. The essential assumption of their derivation is that the 
ttN scattering can be described from the isobar model where the pion is scattered from an 
isobar system which can decay into ttN . If we identify their isobar as of the simple mass 
operator defined by Eq. fl38p . their equation can be schematically cast into the following 
from 

tik',k,E) = v^"P'''\k',k,E) + J dk"v^°P'''\k',k",E)gAAY{k",E)t{k",k,E) (49) 
where the propagator is 

= E-Er,M-EAk)-^%{k,E)+^e ' ^^^^ 

Here ENQ{k) = [{m%)'^ + A;^]^/^ and T.%{k,E) is determined by properly boosting the ex- 
pression Eq. (13T]) . The driving term of Eq. (HOll is determined by the one-nucleon-exchange 
mechanism 

v^"P'^\k',k,E) = [Tmz'J'] [ro(A;')4^'] • (51) 

^ E-EN{k + k>)-E^{k)-E^ik') + ie^ ; o j 

Note that E) and v^°^^^^ in the above equation are defined by the same vertex function 

ro(fc). In the AAY approach, this consistent treatment of the propagator gAAvik^E) and 
the interaction v^"^^^'^ is the consequence of requiring that the scattering equation satisfies 
the TiTiN unitarity condition. In other words, if one replaces the propagator gAAvik^E) 
by g{k,E) of Eq. ( fTSl) . the resulting amplitude from solving Eq. (l49l) will not satisfy the 
unitarity condition. With some derivations, one can also see that the solution of Eq. (H9l) will 
not have a pole at E = rrii^i. The AYY approach simply was not developed to reproduce the 
same analytic structure of the dispersion relations in the unphysical region E < + m^v. 

To further explore the differences with the approaches based on dispersion relations, we 
note that the models considered in this section as well as in the next one solve integral 
equations, such as Eq. f|T7j) . and require form factors to regularize the matrix elements of 
the potential v and the vertex interaction Tq. These form factors can give poles to the 
on-shell scattering amplitudes in the unphysical region of E < rriN + TTt^^ where the nucleon 
pole is identified. More explicitly, if a dipole form is used to parameterize Fq, the on-shell 
matrix element of Eq. (I5U]) becomes 



t{ko,ko,E) 



A2 



k^n + A2 



1 



A^ 



Pa + A' 



(52) 



E-mo-T.%{E) 

where ko is defined hj E = Eisfi^ko) + Ej^i^ko). Thus this amplitude can have pole in the 
region where we define the nucleon pole if A < m^r + m^. This illustrates that the analytic 
structure of the dynamical models deduced from relativistic quantum field theory can not be 
completely consistent with that defined by the dispersion relations of the S-matrix theory. 
In fact, there is no compelling reason to require that they have the same analytic structure. 
In the very extensive literature, as thoroughly reviewed in the textbook of Goldberger and 
Watson [3], the widely used fixed-t dispersion relations in analyzing vrA^ scattering can 
not be derived from relativistic quantum field theory exactly. Historically, the S-matrix 
theory is considered as an alternative to relativistic quantum field theory to study strong 
interactions. There is no rigorous theoretical argument to favor one of them in developing 
phenomenological models to analyze the data. 
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III. METHOD OF UNITARY TRANSFORMATION 



The method of unitary transformation was essentially based on the same idea of the Foldy- 
Wouthuysen transformation developed in the study of electromagnetic interactions. Instead 
of considering the original Lagrangians with bare masses and bare vertex interactions, we 
simply ask how the strong interactions can be described with a phenomenological Lagrangian 
defined by the physical masses and physical coupling constants. It is understood that the 
application of such a phenomenological Lagrangian to calculate any amplitude should drop 
loops associated with one-particle states and vertices which are already absorbed in the 
definitions of physical masses and coupling constants. It is an non-trivial problem to justify 
these rules within the exact theory. But such rules are valid in practice since we will only 
consider leading order terms of a perturbative expansion which will be specified later. This 
means that we assume that we have already solved the one-particle problem within a model, 
such as that defined by Eqs. and this problem will not be dealt with in developing 

reaction models. This is the main difference between the models based on three-dimensional 
reductions described in section II and the model based on the unitary transformation. The 
advantage of the latter is that the unitarity condition can be satisfied trivially; in particular 
in handling the multi-channel multi-resonance reactions. Of course, the price we pay is 
that the connection to the theory of nucleon structure is perhaps more remote than the 
approaches based on three dimensional reductions. 

To illustrate the method of unitary transformation, we again consider the simplest phe- 
nomenological Lagrangian density 

L{x) = Lo{x) + Lj{x) (53) 

where Lq{x) is the usual free Lagrangians with physical masses m^v for the nucleon field ipjsi 
and for the pion field 0^, and 

Lj{x) = 'ipN{x)TN,nN^N{x)(l)Ux), . (54) 

Here TN,nN denotes the physical ttNN coupling (~ fnNN)- It is not the bare coupling Fq 
in Eq. ([1]). The Hamiltonian density H{x) can be derived from Eqs. (!53|) - (!54|) by using the 
standard method of canonical quantization. We then define the Hamiltonian as 

H = I H{x,t = 0)dx. (55) 

The resulting Hamiltonian can be written as 

H = Ho + Hi (56) 

with 

Ho = j dk[EN{k)h\h-^ + E^{k)a\a^] (57) 
Hi = Tn^ttN 

= y" dfcidfc2'^^5(^ - ^1 - ^i)[(rAr,^Ar(A;i - A;2)6t6g^agJ + (c.c.)] (58) 

where 6^ and {b and a) are the creation (annihilation) operators for the nucleon and the 
pion, respectively. For simplicity, we drop the terms involving the anti-nucleon operator d"*" 
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and d. Note that H along with the other constructed generators P, K, and J define the 
instant-form relativistic quantum mechanical description of vrA^ scattering. We will work in 
the center of mass frame and hence the forms of these other generators of Lorentz group are 
not relevant in the following derivations. 

The essence of the unitary transformation method is to extract an effective Hamiltonian 
in a 'few-body' space defined by an unitary operator U, such that the resulting scattering 
equations can be solved in practice. Instead of the original equation of motion H\a) = Ea\a), 
we consider 

H'\a) = E^\a) (59) 

where 

H' = UHU^ 

\a) = U\a) (60) 

In the approach of Sato, Kobayashi and Ohtsubo [3i] (SKO), the first step is to decompose 
the interaction Hamiltonian Hj, Eq. fl58|) . into two parts 

Hi = Hf + Hf (61) 

where Hf defines the process a ^ be with > mf, + rric which can take place in the free 
space, and H^ defines the virtual process with < nib + ^c- For the simple interaction 
Hamiltonian, Eq. fISS]) . it is clear that Hf = and Hf = Hj. 

The essence of the SKO method is to eliminate the virtual processes from transformed 
Hamiltonian H' by choosing an appropriate unitary transformations U. This can be done 
systematically by using a perturbative expansion of U in powers of coupling constants. 
As a result the effects of 'virtual processes' are included in the effective operators in the 
transformed Hamiltonian. 

Defining U = exp{—iS) by a hermitian operator S and expanding U = 1 — iS + ... , the 
transformed Hamiltonian can be written as 

H' = UHU+ 

= U{Ho + Hf + Hf)U+ 

= Ho + Hf + Hf + [Ho,tS] + [Hj,tS] + ^^[[Ho,tS],tS]+--- . (62) 

To eliminate from Eq. ( l62l) the virtual processes which are of first-order in the coupling 
constant, the SKO method imposes the condition that 

Hf + [Ho,iS]=0. (63) 

Since Hq is a diagonal operator in Fock-space, Eq. (!63l) implies that iS must have the same 
operator structure of Hf and is of first order in the coupling constant. By using Eq. (I63l) . 
Eq. fl62|) can be written as 



H' = Ho + H'j, (64) 

with 

H'j = Hf + [Hf, iS] + l [Hf, iS ] + higher order terms . (65) 



10 



Since Hf, Hf, and S are all of the first order in the coupling constant, all processes included 
in the second and third terms of the Hj are of the second order in coupling constants. 

We now turn to illustrating how the constructed Hj of Eq. fl65l) can be used to describe the 
7tN scattering if the higher order terms are dropped. We consider the simple Hamiltonian 
defined by Eqs. fl56l)- fl58l) which gives H[ = and = Tjsf^T,:^. Our first task is to find S 
by solving Eq. (jHS]) within the Fock space spanned by the eigenstates of Hq 

Ho\N)=m^\N) (66) 
Ho\k,p) = iE^ik) + EMip))\k,p} (67) 
Ho\h,k2,p) = {{E^{ki) + E^{k2) + Er,{p))\h,k2,P) (68) 



For two eigenstates / and i of Hq, the solution of Eq. fl63l) clearly is 

{mS)\^) = ^^^^ . (69) 
Ef - Ei 

For the considered Hf = Ftv^ttAt we thus get the following non- vanishing matrix elements 

mi'sm = T.M'^) ^^(,)^;;(^),,^^ Mk + p) (70) 



and 



(72) 



(I ko vr\(tS)\k'v') - T* (k.) -6{k'-k2)6{p'~ki-pj) 

E^{ki) + E^{k2) + En{pi) - E^{k') - En{p') 

" ^^-^^^^^i?.(fcO + i?;v(p.)-i?^(p') ^^^^ 

{ms)\hrk.,P,) = T.Mh)- -5ik-kMv'-k2-m) 



E^{k) + E^ip) - E^{k^) - E^{k2) - En{Pi) 

= TN,nN{k2) „ , . ^rrr^ ^ / n (74) 

En{p) - E^{k2) - En{Pi) 

With the above matrix elements and recalling that Hf = and Hf = Ftv^ttAt for the 
considered simple case, the matrix element of the effective Hamiltonian Eq. ( l65l) in the 
center of mass frame {p = —k and p' = —k') is 

{k\H'j\k') = ^5][((^|r^^^^|J)(/|(z^)|^')-(^l(^'5)l^)(^|riv^-iv|^')] (75) 



The only possible intermediate states are |/) = \N) + \Ti{ki)T[{k2)N{Pi)) . By using Eqs. (I70|) - 
([71]) we then obtain 

{k\H'j\k') = v^'\k,k')+v^''\k,k') (76) 
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where 



1 



E^{k) + ENik) - ruN E^{k') + EN^k') - tun 
1 



n,.N{k'm 



+ 



_EN{k) - E^{k') - E^ik + k') 
1 



EN{k')-E^{k)-EN{k + k') 



^N,7TNik) 



(78) 



Note that v^^'' of Eq. fl77j) is due to the intermediate 'physical' nucleon state state |/) = |A^). 
Here we see an important difference between v^^^ and u^p"'*^) of Eq. fl?I]) for the nucleon-pole 
term which is due to a bare nucleon state within the 3dBS models. There is no bare mass 
m% and energy-dependence in v^'^\ This is a consequence of the unitary transformation 
which eliminates the 'virtual' ttN ^ N process. 

With the above derivations, the effective Hamiltonian Eq. (165|) can be explicitly written 



as 



where 



Ho 
V 



H' = Ho + V 



dk [EN{k)b% + E^{k)aia 



dkdk' [v'^'\k, k') + t;(")(A?, k')]a\hKaj:,h_j:, 



(79) 

(80) 
(81) 



To see further the difference between the models from the unitary transformation method 
and the 3dBS method, let us first recall how the bound states and resonances are defined 
in a Hamiltonian formulation. In operator form the reaction amplitude for a Hamiltonian 
Eq. ([79D is defined by 



t{E) = V + V 



1 



or 



t{E) = V + V 



E-Ho + ie 
1 



-m 



V . 



(82) 



(83) 



E-H' + it 

The analytic structure of scattering amplitude can be most transparently seen by using the 
spectral expansion of the Low equation fl83l) 



{k'\m\k) = {k'\v\k)+Y: 



(fcly|$,)(1'Jl^|A;) 
E-e, 

r(+)\/,Tr(+)| 



r^j,> {k'\v\^'^>'){^'^>'\v\k) 



where Eth is the threshold of the reaction channels, $£. and '^^^) are the discrete bound 
states and the scattering states, respectively. They form a complete set and satisfy 

^l-^"..) = (85) 
H'm^^)) = E'\^^^,^) (86) 



(84) 
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Of course bound state energies are below the production threshold Eth- We now note that 
due to the two-body nature of V defined by Eq. ( IHTl) . Eq. (185|) has the one-nucleon solution 
H'\N) = Ho\N) = ttinIN). But it does not contribute to the second term of Eq. (18^ 
because {7iN\V\N) = 0. Thus the amplitude Eq. does not have a nucleon pole which 
corresponds to a bound state with the mass of the physical nucleon and is formed by the 
physical N and vr of the staring Lagrangian Eq. fl53p . This is consistent with the experiment. 

Here we note that in the Hamiltonian formulation, the amplitude {k'\t{E)\k) depends on 
three independent variables: energy E and momenta k and k'. We can analytically continue 
this amplitude to complex i?-plane for any k and k'. In the complex ii^-plane, the bound 
state poles and unitarity cuts are on the real axis of the physical sheet and the resonance 
poles are on the unphysical sheet. These analytic properties with respect to the energy 
variable E are independent of the momentum variables k and k'. As mentioned at the end 
of section II, the on-shell matrix element (/co|t(£')|/co) could have poles from the form factors 
which are needed to regularize the potential V. Thus the analytic structure of {ko\t{E)\kQ) 
can be different from that of the approaches based on dispersion relations. 

To further see the differences with the 3dBS models and the approaches based on disper- 
sion relations, let us solve Eq. (l82l) by considering only V = v^^\ The matrix element of the 
scattering equation defined by Eq. (!82|) is identical to Eqs. (fT7|) - (fT8|) . With the separable 
form Eq. (!77|) of = v^^\ the above equations can be solved explicitly. The solution is 



T{k,k\E) 



N{k,k',E) 
D{E) 



(87) 



with 



D{E) 



D,{E)Y - Do{E)D2iE) 



(88) 



N{k,k',E) = T*{k') 



Do{E) 



+ 



4 {E{k) - m){E{k') - tun) 
l-D,{E))f 1 1 



+ 



2 \E{k)-m E{k')-mN, 
where E{k) = E^^lk) + Ei\f{k) and 

1 



+ D2{E) 



r7v,.7v(fc) (89) 



DoiE) 
Di{E) 
D2{E) 



q^dq |rjv,^^(g)p- ... 

E — E[q) + le 

(fdq \^ n,-wn{<i)\^^77^ TP( \ , ■ \(TP( \ \ 

2[E — E[q) + ie)[E[q) — m) 

q^dq |rjv^Ar(g)p-— — ^q~r~~T7Tu~\ 

A[E — E[q) + ie)[E[q) — my 



(90) 
(91) 
(92) 



With some inspection, one can see that the amplitude Eq. flHTl) does not have the nucleon 
pole condition D{E = 171^) = for any k and k'. This is what one expects from the spectral 
expansion Eq. fl84|) . If we take the on-shell matrix element E = E{ko) = E{k) = E[k'), one 
then finds 



t{ko,ko,E) = r;^7v(A;o 



dijE) ^ dojE) 
E-m A{E-rn?Y 



+ d2{E) 



tNN\.I^O 



(93) 
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where where 



d«(E) = ^, (94) 
ME) = ^5^, (95) 

The first term in the right-hand side of Eq. fp51) does have a pole at £" = M^r of the dispersion 
relations. But it has additional double poles from the second term as well as from the nNN 
form factor which is often parameterized as a dipole form r,rArAf(^) = (A^/(A^ + A;^))^. It 
should be noted that a pole corresponding to a bound state in a Hamiltonian formulation 
must be for arbitrary k and k'. Thus the pole only from the on-shell matrix element Eq. (p3|) 
is not a ttN bound state with mass ttln- Eq. (p3i) shows again that a dynamical model 
deduced from relativistic quantum field theory does not have, and is not required to have, 
the same analytic structure of the amplitudes in the approaches based on dispersion relations 
in the S-matrix theory. 

To end this section, let us mention that the unitarity condition only requires that an 
acceptable model must have a unitarity cut in the physical region E > + m^v. This 
is trivially satisfied in the model defined by the effective Hamiltonian Eqs. fl80l) - fl8Tl) since 
the interaction V is energy independent. This is an important advantage of applying the 
method of unitary transformation to develop a multi-channels multi-resonances reaction 
models for investigating meson-nucleon reactions in the nucleon resonance region, as devel- 
oped in Ref . j3| ■ In a model with an energy-dependent V such as the AAY model the unitarity 
condition is much more difficult to satisfy, and the analytic continuation of the scattering 
t-matrix defined by Eqs. flMl) to complex ii^-plane is in general much more involved. 



IV. TIME-ORDERED PERTURBATION THEORY 

Treating the Hamiltonian Eq. fISBl) in time-ordered perturbation theory [3], the matrix 
elements of the transition operator can be represented by a series expansion defined by all 
diagrams containing an incoming and outgoing pion-nucleon state 

{nN\t{E)\nN) = {nN\Hi- -Hi\nN) 

E — Hq + ie 

+ (nNlHi Hi Hj ^ HAnN) 

^ ' E-Ho + ie E-Ho + ie E-H^ + ie ' ' 

+■■■. (97) 

In the approach of the Julich group js], 0], the partial sum of this series is written as a 
three-dimensional integral equation which takes the same form of Eqs. fl 171) - (1181) . In the 
simple model defined by Eq. fl57j) - fl58l) . the resulting potential can be schematically written 
as (dropping anti-nucleon terms) 

v{k, k', E) = v^'^k, k', E) + ^("^A?, k\ E) (98) 
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where 



E 



m 



N 



E - EN{k + k') - E^{k) - E^{k') 



N,ttN 



(k) 



(99) 

(100) 



As can be readily seen the term v^{E) has the singularity of the vrvrA^ cut, and they depart 
from the starting Hamiltonian Eq. (!57 l) -( !58l) by using the bare mass to define the s- 
channel term v^^\ The above potential has the same form of the matrix element of that 
defined by Eqs. (l20 l) -( l2Ti) . Thus the interpretation of their analysis of nucleon pole term is 
similar to what described in section II. 



V. SUMMARY 

In this paper, we have examined three methods for constructing meson-nucleon reaction 
models from relativistic quantum field theory. For the models based on the three-dimensional 
reductions of Bethe-Salpeter equation and the time-ordered perturbation theory, the driving 
terms of the resulting three-dimensional scattering equations in general contain a nucleon 
pole term determined by a bare nucleon Nq. We show that the commonly used procedure 
of imposing the nucleon pole condition to fix the bare nucleon parameters is related to the 
assumption that the nucleon is a bound state made of a bare core A'^o and meson cloud. To 
correctly implement this nucleon substructure into the scattering equation, it is necessary 
to consider TrvriV unitarity condition as achieved within the model of Aaron, Amado and 
Young (l2|. 

We have given a pedagogical and explicit explanation of the method of unitary trans- 
formation which has been applied in recent years to investigate meson-nucleon reac- 
tions 0, 3, m, 16, 17|. Since only physical nucleons and pions are the basic degrees of 



freedom of the derived effective Hamiltonian, the resulting vrA^ amplitude does not have a 
nucleon pole aX E = Mn in the complex i?-plane. This is due to the fact that the one- 
nucleon problem is decoupled from the two-particle problem by the unitary transformation 
and the resulting vrA^ potential is energy independent. There is no nN interaction due to the 
propagation of a hare nucleon in this formulation. We explain how this can be understood 
from the general principles of a Hamiltonian formulation of reactions. 

We also show that the scattering amplitudes from these three methods do not have 
the same analytic structure of the amplitude from the approach based on the dispersion 
relations of the S-matrix theory. Even one imposes the nucleon pole condition to make the 
connection to the S-matrix theory, the constructed models can have poles from the form 
factors which are needed to regularize the potentials for solving the resulting scattering 
equations. We emphasize that there is no compelling and rigorous theoretical argument to 
request that the constructed models should have the analytic structure of the dispersion 
relations. There exits no rigorous derivation of the dispersion relations of S-matrix theory 
from relativistic quantum field theory 0]. Historically, the S-matrix theory is considered as 
an alternative to relativistic quantum field theory to study strong interactions. Either one 
of them is a good starting point for developing phenomenological models for analyzing the 
data. For investigating multi-channels and multi-resonances reactions, the models deduced 
from relativistic field theory appear to be more practical. 
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